Role of Impurities in Core Contribution to Vortex Friction 
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We present a microscopic calculation of core states contribution to vortex friction in a type II 
superconductor at zero temperature. In weak impurity limit, a perturbative calculation leads to a 
vortex friction proportional to the normal state resistivity. In dirty limit, vortex friction due to core 
states is determined only by the Fermi energy and the energy gap, similar to the phenomenological 
results of Bardeen and Stephen, but disagrees with the results obtained under the relaxation time 
approximation. In addition, an explicit demonstration is given on the the insensitivity of transverse 
force to impurities. 
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The importance of vortex dynamics in a superconduc- 
tor has long been realized and been under continuous 
and intensive studies S. In the early work by Bardeen 
and Stephen, it was found that in dirty limit and at zero 
temperature, vortex friction is determined by the Fermi 
energy and the energy gap of the superconductor. [Q This 
phenomenological result has ample experimental support 
||. Subsequent microscopic derivations, however, have 
reached a conclusion that the vortex friction should ap- 
proach zero in dirty limit. Meanwhile, these calculations 
also give vanishing a transverse force on a moving vortex 
in the same limit. B 

The magnitude of transverse force on a moving vor- 
tex has attracted much attention in recent years. The 
Berry phase acquired by a moving vortex shows that the 
Mangus force is the only transverse on the vortex at zero 
temperature irrespective whether or not impurities are 
presented. || This observation has been confirmed by a 
direct measurement of the transverse force in dirty su- 
perconductors ||. It has led to a re-examination of the 
approximations used in the earlier microscopic deriva- 
tions of vortex dynamics. It has been found that the 
relaxation time approximation always leads to incorrect 
results when used in force calculations by either force cor- 
relation or force balance conditions. Q| The calculations 
which give a vanishing vortex friction in dirty limit, in 
disagreement with the Bardeen-Stephen result, are ex- 
actly those using relaxation time approximations. 

The problem now becomes how to calculate vortex fric- 
tion from a microscopic theory in presence of impurities 
without the improper use of the relaxation time approx- 
imation. Solving this problem is the main purpose of 
our paper. In the present paper we will report a set 
of results based on the BCS theory which shows that 
core states contribution to vortex friction arises and sat- 
urates as the superconductor evolves from weak impurity 
to dirty limit, determined by the energy spectral property 
of the BCS Hamiltonian. An explicit demonstration of 
the insensitivity of the transverse force to impurities will 
also be given here, providing another way to concretely 
realize the results from the Berry phase calculation. 

We consider an isolated rectilinear vortex at x v of 
length L in a BCS superconductor. For our purpose of 



calculating core states contribution to vortex dynamics, 
we will not include the vector potential of a magnetic field 
explicitly because its effect on the core state structure is 
negligible. The vortex velocity is given by v„ = x v . The 
distribution of impurities is assumed to be homogeneous 
and the variation of impurity potential only appreciable 
on a scale much smaller than size of the vortex core. If 
the vortex is allowed to move slowly against the ionic 
lattice background, to the leading order of the vortex 
velocity, according to the Ehrenfest theorem, the cor- 
responding force acting on the vortex has the form j|] 
F = Bv v x z — r/v v . Within the BCS theory, the trans- 
verse coefficient B of the transverse force is explicitly ex- 
pressed in terms of the transverse force-force correlation 
function MM 



B = 
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and the vortex friction coefficient r\ in terms of the lon- 
gitudinal force-force correlation function 
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Here f k = 1/(1 



) is the Fermi distribution func- 



tion. The wavefunctions {^(x)} and the corresponding 
eigenvalues {Ek} are determined by the Bogoliubov-de 
Gennes equation, TL ^k(x) = Ek ^k(x) , with *ff k (x) — 
u k (x) 
v k (x) 

11 A ' " " (h 2 /2m)V 2 - fi F + V{x), 



, and the system Hamiltonian given by TL 



Here H 



A* -H* 

and V(x) the impurity potential. The order parameter 
is determined self-consistently 

A(a:) = g < YP- ~ Vk)u k {x)v* k {x) > , 



where the impurity average < ... > is implied. The am- 
plitude of the order parameter |A| is the energy gap, 
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and g the strength of attractive interaction between elec- 
trons. Since the impurity distribution is homogeneous, 
its average leaves the vortex position the only reference 
point in the system. We expect A(x, x v ) = A(x — x v ) = 
|A(x — x v )\e t8 ^ x ~ Xv > with 9 the singular phase. This gives 
a convenient specification of the vortex position in the 
calculation, which we will adopt in this paper. 

We will first evaluate the transverse force to show its 
insensitivity to impurities. We demonstrate that B can 
be expressed in a manner without explicit dependence on 
the system Hamiltonian. Because TL is Hermitian, all its 
eigenstates form a complete and orthonormal set, that 
is, J d 3 x^l(x)^ k ,(x) = 5 k ,k> , and J2 k ^k(x)^l(x') = 
1 , with $! k (x) = (u^.(x),vl(x)) and the wavefunction is 
normalized to 1 over a cylinder of radium R and length 
L. The thermodynamic limit requires that R — > oc. For 
states k ^ k' we have 



d 3 x yl(x)(V v H)V k ,{x)/(E k , - E k ) 

d 3 x y{(x)V v if kl {x) = - [ d 3 x V v y{(x)y k ,(x) , (3) 



where a detailed solution of these states from the 
Bogoliubov-de Gennes equation with a vortex can be 
found in Ref. pl| . Eq.(3) is obtained by taking gradi- 
ent V„ with respect to H\^ k >) = E k >\^ k >) and {^ k \TL = 
E k {^ k \, then multiplying from left or right by {^ k \ or 
l^fc') respectively. With the aid of Eq.(3) and the com- 
pleteness of the eigenfunctions, the transverse coefficient 
can be expressed as a summation over individual state 
contributions, 



B = ih^2 / d 3 x{f k V v ii k (x) x V v u* k (x) 
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Now we are ready to explicitly demonstrate that the coef- 
ficient B for the transverse force is insensitive to impuri- 
ties as first observed through the Berry phase calculation. 

In a clean superconductor we can use the replacement 
of V v — > —V in Eq.(4). Using this replacement and the 
definition of number current, 



- ^ {fkU* k Vu k + (1 - f k )v k Vv* k } + c.c. , 



2 ^ 

k 



Eq.(4) becomes 



B = h J d 3 x z ■ (V x j) 
= Lh <z> cfl • j 
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At zero temperature, f k = 1(0) for E k < 0(> 0), the cal- 
culation of p s is straightforward. It is equal to the total 



tion: HH^(x) = E , *(x) 



then 



superfluid number density po — J2 k E k >o \ v k(\ x — x v\ — ► 
oo)| 2 , number of Cooper pairs per unit area. Here we 
have used the relations between the positive and neg- 
ative energy states of the Bogoliubov-de Gennes equa- 

v*(x) 
-u*(x) 
TL $>(x) = -E . 

In the presence of impurities, the replacement V„ — > 
—V cannot be directly used in Eq.(4) because of the 
implicit impurity dependence. Instead, we expand ^ k 
in terms of eigenfunctions {<!>;} of TL, a corresponding 
Hamiltonian to TL without the impurity potential V(x) 
and with impurity averaged A, 
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Here {a kll } and {5 kll } are the modulus and phases of the 
expansion coefficients. They are functions of the vortex 
position as well as the positions of impurities. We assume 
these coefficients to be described separately by two inde- 
pendent random matrices, making use of the randomness 
of impurities. |f2] ] TL still has a dependence on impurities 
through the impurity averaged order parameter A: 



n 



H 

A* 



A 

-H* 



(7) 



with H = -(h 2 /2m)V 2 - p F , TL<$>i = Ei^ t and $/ = 
_ j . Since form a complete set, the expan- 

sion coefficients {a k ie lSkl } form a unitary matrix, and 

E 



l or k u kl 



= 1. For <i>; we can use the replacement 
V„ — > —V. We remind here that although away from 
the vortex core the value of the order parameter is the 
same as that in the clean case, guaranteed by the An- 
derson theorem, in the core regime the profile of A is 
different, corresponding to a different coherence length. 
The core size will be smaller with impurities presented. 
Using Eq.(6), Eq.(4) becomes 

B = -ih < [ d 3 x{f k V v (a k ie lSk 'ui{x)) x 

k 1,1' J 

V, {a kl ,e- i5 ^u* v {x)) - (1 - / fc )V„ (a kl e iS *%(x)) x 
V v (a kl >e- iS »i'vp(x))}-z> . (8) 

Here < ... > stands for the impurity average over the 
expansion coefficients. Then, 

B = -i^ y^y^ < a 2 kl > / d 3 x {f k V v ui(x) x V v u*(x) 
t. i J 
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Terms containing derivative to vortex position in the ex- 
pansion coefficients have been averaged to zero. Now the 
replacement of V u — > — V can be used, and following the 
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procedure leading to Eq.(5) we have, at zero tempera- 
ture, the desired result B = 2irhpQ , because 

< a ki > \vi(\x-x v \ -> oo)| 2 

k,E k >0 I 

= 2J < \ v k(\ x - x v\ -» oo)| 2 >= Po • (10) 

The above second equality is the Anderson theorem. 
Eq.(10) can also be reached from the envelop wavefunc- 
tion argumentation fl3| . 

In the following we turn to the calculation of the core 
states contribution to vortex friction at zero tempera- 
ture. We start with a clean superconductor at zero tem- 
perature. In this limit core state energies are completely 
labeled monotonically by a half integer fi associated with 
the azimuthal angle in the eigenfunction jfljand the 
core energy spectrum is discrete. In the transition matrix 
element shown in Eq.(2), only those between neighbor- 
ing fi are non-zero. In addition, because of the factor 
(fk — fk') in Eq.(2), only the two states closest to the 
Fermi surface, one above and one below, need to be con- 
sidered. However due to the finite energy level spacing, 
their contribution to the vortex friction is zero, because 
of the delta function in energy in Eq.(2). Therefore, at 
zero temperature and in the clean limit there is no core 
states contribution to vortex friction. 

In the presence of a finite impurity potential V(x), the 
energy spectrum of core states no longer has a monotonic 
dependence on [i. Transition elements in Eq.(2) are no 
longer confined to neighboring fi states. The impurity po- 
tential also brings eigenenergies of some of the core states 
close to each other so that it generates a quasi- continuous 
energy spectrum after impurity average, necessary for a 
finite vortex friction according to Eq.(2). In the weak 
impurity potential limit, their contribution to the vortex 
friction can be calculated perturbatively. We will show 
that vortex friction due to core states is proportional to 
the normal state resistivity. 

To explicitly exploit the perturbation method, it is 
more convenient to use V^> k (x) instead of V v ^f k (x) in 
Eq.(2). It may be done in the following manner. The 
system Hamiltonian may be written as Ji = H + W, 
with TL defined by Eq.(7), and 



w = 



V(x) 
-V{x) 



Next we use 



(VH) = VH+VH' = -(V v n)+VW = ~{V V H) + \7H' , 

by noting that Ti. and V v Ti. have same dependence on x v , 
in addition to 
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and Eq.(3) to obtain 
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The effect of impurities is now described by VW. With 
Eq.(ll), the transition element is given by 
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For the leading order of contribution, we use the unper- 
turbed and A. The first term and the cross terms 
in Eq.(12) will not give any contribution to the vortex 
friction because of the discreteness of {Ey — E k ) and 
the delta function factor 5(E k — Ey) in the vortex fric- 
tion coefficient rj. The contribution to 77 comes from the 
last term. Without losing generality, we assume that the 
impurity potential has a length scale smaller than the 
coherence length such that we may effectively describe it 
by a delta potential V(x) — J2i VoS 3 (x — X4), with Vq car- 
rying a random sign and xi randomly distributed. Using 
the wavefunction obtained in Ref. JllJ] , we obtain the last 
term of Eq.(12), 



< 



d 3 x^\{x)VH'^ k >{x) 
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Here rii is the impurity concentration. Hence the vortex 
friction coefficient is 



7T / kpL 

n = —fi 
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Here eo = A^/Ep is the core level spacing, ^kp — 
Ep/Aoo, Ep = (hkp) 2 /2m, and Aoo is the energy gap 
far away from the vortex core. In Eq.(13) kpL/(2Treo) = 
n c (E) is the approximated number of core states per unit 
energy near the Fermi surface, arose from the summation 



J2S(Ek-E k ,) 



fk - fk' 
E k i — E k 



dE k dE k , 5{E k -E k ,) 



fk - fk' 

E k > — E k 
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We now make connections to the normal states transport 
coefficients. Note that for the normal state the electronic 
relaxation time and the electron scattering cross section 
have the following relations, t^ 1 = niVpUtr , with the 
transport cross section cr tr = J dfl (1 — cos 8)\V (9)\ 2 . 
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@ Here V{6) = -~m/2irh 2 J d 3 xV (x)e~ l ^ r , with q = 
k — k', 8 is the angle between k and k', and vp = hkp/m. 
For our choice of potential, r^} = niVp(m/2iTTi 2 ) 2 VQ 4n . 
Expressed in r tr for the transition element, the vortex 
friction coefficient becomes 



3 mn e (irt; 2 L) 



(14) 



with the electron density n e — k F /(3ir 2 ) . The scatter- 
ing time r tr may be related to the normal state residual 
resistivity by p = raj ' (n e e 2 T tr ), and it can be measured 
independently. Eq.(14) shows that the vortex friction in 
the weak impurity limit increases with the normal state 
resistivity We show below that this increase eventually 
saturates in the dirty limit. 

Using Eq.(6), we expand localized states in terms 
of {$;}, the eigenfunctions of H. Because \7 V H = V^H, 
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In the above derivation, we have ignored the localiza- 
tion effect which depresses the density of state, or the su- 
pcrfiuid density. We justify our assumption here. There 
are three energy scales involved in the derivation of vor- 
tex dynamics, the Fermi energy Ep, the energy gap A^, 
and the core level spacing A^/Ep. The effect of im- 
purities on vortex dynamics is believed appreciable at 
TtrA^/hEF < 1 {§, and the equal of above Eq.(14) to 
(16) suggests that the impurity starts to be effective at 
Tt r A 00 /h(A 00 /Ep) 2 ~ 1. They indicate that the impu- 
rity effect on vortex friction occurs at a rather weak level, 
determined by the smallest energy scale in the problem. 
The dirty limit is given by A^/Ep < T^A^/h < 1. 
The localization effect is only pronounce in the extremely 
dirty limit, the localization regime, when r tr Ep/h < 1. 
|l5| Because Aoo /Ep « 1, away from the localization 
regime the depression of density is indeed negligible. The 
un-depressed electronic density applies. 

In conclusion, we have found that the impurity aver- 
age smoothes the core state energy spectrum, a process 
necessary for vortex friction. We have obtained that the 
core states contribution to vortex friction is proportional 
to the normal state resistivity in weak impurity limit and 
saturates in dirty limit. 
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Here t c — kp/2. In dirty limit, all core states of TL are 
combined to make the core states of 7i, that is, there is a 
uniform distribution for aa in terms of indices of both I 
and k. The normalization condition for the wavefunction 
requires then < >— 1/N C , with the total number 
of core levels N c w 2(A oc /eo)(kpL/2ir), same as in clean 
superconductors. This is implied by the Anderson theo- 
rem that there is no change in the number of extended 
states as the superconductor evolves from clean to dirty 
limit. Hence we can treat the core states and extended 
states separately, and ignore the mixing between those 
two groups of states. The average value of the transition 
element is now 
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d 3 x 9Ux) (V,H)f fc -W 



>= 



4e 2 \t c \ 2 
N c 



and the vortex friction coefficient is then 

tt fk F L\ 2 4e 2 \t c \ 2 37T% Aoo 

T) = —n — ^ = nn e L 

' 2 \2ireaJ N c 8 Ep 
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This result is similar to the value obtained in Ref. [jjj. 
Hence its microscopic justification has been provided for 
the first time. The magnitude of this vortex friction con- 
tribution is smaller than the transverse force by a factor 
of Aoo/Ep. 
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